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Abstract 

Wc have measured the branching ratio BR((/) 7r'^7r'^7) with the KLOE detector 
using a sample of '^5 x 10'' (p decays. mesons are produced at DA$NE, the 
Prascati ^-factory. We find BR{(l)^Tr'^Tr^-f)= (1.09 ± O.OS^tat ± O.OSsyst) x IQ-^. We 
fit the two-pion mass spectrum to models to disentangle contributions from various 
sources. 

Key words: e+e" collisions, (p radiative decays, scalar mesons 
PACS: 13.65.+i, 14.40.-n 



The decay — > 7r°7r°7 was first observed in 1998 [1]. Only two experiments 
have measured its rate [2,3]. The measured rate is too large if >/o(980)7, 
with /o TT^TT^, were the dominating contribution and /o(980) is interpreted 
as a qq scalar state [4,5]. Possible explanations for the /o are: ordinary qq 
meson, qqqq state, KK molecule [6-8,4]. Similar considerations apply also to 
the ao(980) meson. The decay — > n^n^^ can clarify this situation since both 
the branching ratio and the line shape depend on the structure of the /q. 
We present in the following a study of the decay — > n^iT^'-f performed with 
the KLOE detector [9] at DA$NE [10], an e'^e^ collider which operates at a 
center of mass energy W=M^ ~ 1020 MeV. Data were collected in the year 
2000 for an integrated luminosity L^t ~ 16 pb~^, corresponding to around 
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5 X 10^ 0-meson decays. 



The KLOE detector consists of a large cylindrical drift chamber, DC, sur- 
rounded by a lead-scintillating fiber electromagnetic calorimeter, EMC. A su- 
perconducting coil around the EMC provides a 0.52 T field. The drift cham- 
ber [11], 4 m in diameter and 3.3 m long, has 12,582 all-stereo tungsten sense 
wires and 37,746 aluminum field wires. The chamber shell is made of carbon 
fiber-epoxy composite and the gas used is a 90% helium, 10% isobutane mix- 
ture. These features maximize transparency to photons and reduce — > Kg 
regeneration and multiple scattering. The position resolutions are (7a;j^~150 
/im and (T^~ 2 mm. The momentum resolution is (T{p±)/p± ~ 0.4%. Vertices 
are reconstructed with a spatial resolution of ~3 mm. The calorimeter [12] is 
divided into a barrel and two endcaps, for a total of 88 modules, and covers 
98% of the solid angle. The modules are read out at both ends by photomulti- 
pliers; the readout granularity is ~4.4x4.4 cm^, for a total of 2440 cells. The 
arrival times of particles and the positions in three dimensions of the energy 
deposits are obtained from the signals collected at the two ends. Cells close in 
time and space are grouped into a calorimeter cluster. The cluster energy E 

— * 

is the sum of the cell energies. The cluster time T and position R are energy 



weighted averages. Energy and time resolutions are ge/E — b.7%/ JE (GeV) 



and at — 57 ps/y E (GeV) © 50 ps, respectively. The KLOE trigger [13] uses 
calorimeter and chamber information. For this analysis only the calorimeter 
signals are relevant. Two energy deposits with > 50 MeV for the barrel and 
E > 150 MeV for the endcaps are required. 

Prompt photons are identified as neutral particles with (3 — 1 originated at 
the interaction point requiring \T — R/c\ < min(5(TT, 2 ns), where T is the pho- 
ton flight time and R the path length; ctt includes also the contribution of the 
bunch length jitter. The photon detection efficiency is ~ 90% for E^=20 MeV, 
and reaches 100% above 70 MeV. The sample selected by the timing require- 
ment contains a < 1.8% contamination due to accidental clusters from machine 
background. 

Two amplitudes contribute to — >• TT^ir'^'y: (p^S'j, S^tt^tt^ {Sj) and (p^p^ir^, 
p°^7r°7 (pyr) where is a scalar meson. The event selection criteria of the 
(j) — > TT^TT^^ decays (77777) have been designed to give similar efficiencies for 
both processes. The first step, requiring five prompt photons with E^^ >7 McV 
and 6 > ^^min = 23°, reduces the sample to 124,575 events. The background 
due to 4>^KsKl is removed requiring that £^tot = Z^s E^,i and ptot = Hb'P'y,i 
satisfy ^tot >800 MeV and |ptot| <200 MeV/c. We are left with 15,825 events. 
Other reactions which give rise to background are: e"'"e~— >a;7r'^— >7r°7r°7 (a;7r), 
0— >777r°7— >57 (77777) and 0— ^777— >37r°7 (777) with 2 undetected photons. 

A kinematic fit (Fitl) requiring overall energy and momentum conservation 
improves the energy resolution to 3%. Photons are assigned to 7r°'s by mini- 
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Fig. 2. Efficiency vs vr vr invariant mass for (p vr vr 7 events. 

mizing a test x^-function for both the 71777 ^ind ujtt cases. For the ujh case we 
also require Mj^^ to be consistent with M^. The correct combination is found 
by this procedure 89%, 96% of the time for the 77777, ^'^ case respectively. 
Good agreement is found with the Monte Carlo simulation, MC, for the dis- 
tributions of the ^-iid of the invariant masses. A second fit (Fit2) requires 
the masses of 77 pairs to equal M^r. 

The e+e~ utc'^ 77°77°7 background is reduced rejecting the events satis- 
fying x^/ndf < 3 and AM^^ = |M^^ — M^\ < using Fit2 in the a;77 
hypothesis. Data and MC are in good agreement (Fig. l.a-b). The — >■ 77°77'^7 
events must then satisfy x^/ndf < 3 for Fit2 in the 77777 hypothesis. We also 
require AM^^ = \M^^ — Mj^\ < 5cr^ using the photon momenta of Fitl. The 
efficiency for the identification of the signal is evaluated applying the whole 
analysis chain to a sample of simulated (p^Sj, S^h^tt^ events with a 77^77° 
mass (m) spectrum consistent with the data. We use the symbol M^^r to denote 
the reconstructed value of m. The selection efficiency as a function of M^^r is 
shown in Fig. 2. The average over the whole mass spectrum is e-j^T^^ = 41.6%. A 
similar efficiency function is obtained for the process — >■ p°77° with p° — ^ 77^7. 
Fig. 3 shows various distributions for the 3102 events surviving the selection 
together with MC predictions. The angular distributions prove that is the 
dominant process. The rejection factors and the expected number of events 
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Fig. 3. Data-MC comparison for (p — > 7r'^7r'^7 events after lot: rejection: (a) x^/^idf; 
(b) (M^^ — Mjr)/c77r with x^/ndf < 3; (c, d) angular distributions with all analysis 
cuts applied. 6 is the polar angle of the radiative photon, ■0 is the angle between 
the radiative photon and vr" in the ir^n^ rest frame. 



Table 1 



Background channels for 


<j) — > 7r°7r'^7. 




Process 


Rejection Factor 


Expected events 


e~^e~ — > ujt:^ ir^Tr^j 


8.7 


339 ± 24 


(f) — > r/7r'^7 777r'^7 


4.0 


166 ± 16 


~* ^7 ~^ 7r'^7r'^7r'^7 


5.9 X 10^ 


159 ± 12 



for the background processes are listed in Tab. 1 [14-16]. After subtracting 
the background 2438±61 — > 7r°7r°7 events remain. Their Mj^^^ spectrum is 
shown in Fig. 4. 

The (f) 7r'^7r°7 branching ratio, BR, is obtained normalizing the number of 
events after background subtraction, N — B, to the (p cross section, a{(f)), to 
the selection efficiency and to Ljnt: 

BR(0 vr^A) = — i— (1) 



The luminosity is measured using large angle Bhabha scattering events. The 
measurement of cr(0) is obtained from the ^ 777 — 777 decay in the same 
sample [15]. We obtain: 

BR(0 TtOA) = (1.08 ± O.OSstat ± O.OSsyst ± 0.04norm) X 10"^ (2) 
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Table 2 

Uncertainties on BR(^ — > 7r'^7r'^7). 


Source 


Relative error 


Statistics 


2.5% 


Background 


1.3% 


Event counting 


2.3% 


Normalization 


3.7% 


Total 


5.2% 



The contributions to the uncertainties are listed in Tab. 2. Details can be 
found in Ref. [16]. 

In order to disentagle the contributions of the various processes and to deter- 
mine the normalized differential decay rate, dBR/dm = (l/r)dr/dm, we fit 
the data to a mass spectrum f{m). This spectrum is taken as the sum of ^7, 
pTT and interference term, /(m) = fs'y{'^) + fpn{^) + /mtl'^)- The scalar term 
is [17]: 

fsj{m) = — — ^,1^ — . (3) 



The 4>^Sj process is estimated by means of a K^K loop for the /q: 



where g^K+K- and gf^K+K- are the couplings and g{m) is the loop integral 
function. 

A recent measurement [18] reports the existence of a scalar a with = 
(478l^^±17) MeV and = (324l^^±21) MeV. If we include the contribution 
of this meson, its decay rate is given by [19]: 

^ , , e'gl., 1 f Ml-m^ Y 



where g^^-y is a point-like 0(T7 coupling. 
Finally, F^TrOTrO is related to Tsn+Tr- by: 



1 a J. - 4M^ 

F5.o,o (m) = - F5.+,- (m) = ^fi^ 1/1 f. (6) 

2 Szirm \ 
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Fig. 4. Observed spectrum of invariant mass before (a) and after (b) back- 

ground subtraction. 

Table 3 

Fit results using /o only, Fit (A), and including the a, Fit (B). 





Fit (A) 


Fit (B) 


xVndf 


109.53/34 


43.15/33 


M/o (MeV) 


962 ±4 


973 ± 1 


5?o/.+i.-/(4vr) (GeV2) 


1.29 ±0.14 


2.79 ±0.12 


9%K+K-/9%n+7T- 


3.22 ±0.29 


4.00 ±0.14 


94icr'y 




0.060 ± 0.008 



For the inverse propagator, Ds, we use the formula with finite width correc- 
tions [17] for the /o and a Breit Wigner for the a. The parametrization of 
Ref. [20] has been used for the pvr and the interference term. 

The observed mass spectrum Sobsl^-n-Tr) is fit folding into the theoretical shape 
experimental efficiency and resolution after proper normalization for ct(0) and 
Lint- Two different fits have been performed varying /^^(m): in Fit (A) only the 
/o contribution is considered while in Fit (B) we also include the contribution 
of the a meson. The mass and width of the a were fixed to their central values. 
If the normalization of the pvr term is left free during fitting, its contribution 
and the associated interference terms turn out to be negligibly small. When 
BR(0 -> p°7rO 7r°7r°7) is fixed at 1.8 ■ 10"^ as in Ref. [20], the xVndf 
increases by more than a factor of 2. The fits without the pvr contribution are 
shown superimposed over the raw spectrum in Fig. 4.b. 

The result of the fits are listed in Tab. 3. In Fit (A) we use as free parameters 
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^h- 9%K+K- and the ratio g%K+K- 1 9%n+7.-- The fit gives a large xV^df; 
integrating the theoretical spectrum a value BR(0 fgj Tc^n^j) = (1.11 ± 
0.06stat+syst) X 10"^ is obtained. 

A much better agreement with data is given by Fit (B), where we add as 

a free parameter also the coupling gff,^^. The negative interference between 
the /o and a amplitudes results in the observed decrease of the n^^n^j yield 
below 700 MeV. Integrating over the theoretical a and /o curves we obtain 
BR(0 ^ a7 ^ = (0.28 ± O.O^tat+syst) x IQ-^ and BR(0 ^ /07 ^ 

7r07r°7) = (1.49 ± 0.07stat+syst) X IQ-"^. Multiplying the latter BR by a factor 
of 3 to account for /o — > tt+tt"" decay, the BR(0 — > /07) is determined to be: 

BR(0 ^ /07) = (4.47 ± 0.21,tat+syst) ■ 10-^ (7) 



The values of the coupling constants from Fit (B) are in agreement with 
those reported by the SND and CMD-2 experiments [2,3]. The coupling con- 
stants differ from the WA102 result on /q production in central pp coUisions 
(9%K+K-/9%n+-K- = = 1.63 ± 0.46) [21] and from those obtained 

when the /o is produced in Df — > 7r"'"7r~7r"'' decays [22], where qk is consistent 
with zero. 

In order to allow a detailed comparison with other experiments and theoretical 
models, we have unfolded 5'obs(-^7r7r)- For each reconstructed mass bin, the ra- 
tio between the theoretical and the smeared function, SF(Mt^t^), is calculated. 
The dBR/dm is then given by: 

dBR ^ ,gobs(M^^) 1 

dm SF{Mt^t^) Lint X 



The value of dBR/dm as a function of m is given in Tab. 4 and shown in 
Fig. 5. Integrating over the whole mass range we obtain: 

BR(0 ^ ttOtt^) = (1.09 ± 0.03,tat ± 0.03,yst ± 0.04„„,„,) x 10"^ (9) 



which well compares with the result obtained correcting for the average selec- 
tion efficiency (Eq. 2). If we limit the integration to the /o dominated region, 
above 700 MeV, we get: 

BR(0 ^ 7r°7r%; m > 700 MeV) = 

(0.96 ± 0.02stat ± 0.02syst ± 0.04norm) X 10"^ 

which is in agreement with our previous measurement in the same mass range 
[23]. 
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Fig. 5. dBR/dm as a function of m. Fit (B) is shown as a solid line; individual 
contributions are also shown. 

Table 4 

Differential BR for (j) tt^tt^^. m is expressed in MeV while dBR/dm is in units 
of 10*^McV^^. The errors listed arc the total uncertainties. 



m 


dm, 

dm. 


m 


c/BR 
d.m 


290 


2.0 ±2.9 


670 


11.2 ± 1.9 


310 


2.2 ± 1.4 


690 


11.0 ± 1.9 


330 


3.0 ± 1.5 


710 


12.5 ± 1.9 


350 


0.9 ±1.3 


730 


14.0 ± 2.0 


370 


2.9 ± 1.4 


750 


17.3 ± 2.3 


390 


2.2 ± 1.3 


770 


17.0 ±2.4 


410 


1.4 ± 1.1 


790 


19.4 ±2.5 


430 


1.8 ±1.0 


810 


27.4 ± 3.1 


450 


1.9 ±0.8 


830 


29.2 ±3.2 


470 


1.1 ±0.5 


850 


30.6 ± 3.2 


490 


0.5 ±0.2 


870 


41.7 ±3.8 


510 


0.2 ±0.1 


890 


39.6 ±3.6 


530 


0.3 ±0.2 


910 


44.6 ± 3.8 


550 


1.3 ±0.5 


930 


53.6 ±4.4 


570 


3.3 ± 1.5 


950 


47.2 ±4.3 


590 


2.1 ±3.6 


970 


64.7 ±5.3 


610 


3.7 ±4.7 


990 


22.0 ± 2.5 


630 


4.2 ±3.7 


1010 


0.2 ±0.1 


650 


7.0 ± 1.7 







In a separate paper [14], we present a measurement of BR(0 — > ao7), together 
with a discussion of the implications of /o and qq results. 
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